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(54) Method for determining the attenuation of a PCM signal over a digital channel 



(57) A method of determining digital channel atten- 
uation; comprising the steps of: receiving a known train- 
ing sequence of PCM codes, which PCM codes are 
subjected to the attenuation within the digital channel; 
quantizing the received known training sequence of 
PCM codes according to a predetermined thresholding 
procedure; identifying identical PCM codes created as a 
result of the thresholding procedure; and, determining 
the attenuation of the digital channel based upon the 
identification of identical PCM codes. A method is also 
disclosed for determining a digital channel PCM code 
transformation comprising receiving a known training 
sequence of PCM codes, which PCM codes are sub- 
jected to the PCM code transformation within the digital 
channel, quantizing the received known training 
sequence of PCM codes according to a predetermined 
thresholding procedure, and determining the transfor- 
mation of transmitted codes to those received . A 
method is also disclosed for improved echo cancellation 
in a communications network having an analog and a 
digital modem, comprising saving codes transmitted 
from the digital modem to the analog modem for echo 
cancellation, transforming, by a mapping table, codes 
transmitted from said digital modem to codes received 
by the analog modem, and, using the received codes as 
a reference signal for cancellation of echo. A method of 
improved spectral shaping using a transmit shaping 
transfer function in a communications network having 
an analog and a digital modem, comprising, transform- 
ing, by a mapping table, codes transmitted from the dig- 
ital modem to codes received by the analog modem, 



using the received codes for transformation to their lin- 
ear value equivalent representations, and, applying the 
linear value representations to the transmit shaping 
transfer function. 
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Description 

Field of the Invention 

This invention pertains generally to modem technology, more specifically to PCM modem technology, and more 
specifically, to a method for the discovery of the digital attenuation from a set of received PCM samples corresponding 
to a known set of transmit PCM samples. 

Background of the Invention 

This invention provides a method for the discovery of the digital attenuation from set of received pulse code modu- 
lation (PCM) samples corresponding to a known set of transmit PCM samples. The principal feature set used is the 
knowledge of which received PCM codes have become indistinguishable (i.e., are identical) as a result of the attenua- 
tion mapping. An alternative but equivalent feature set would be use of the absent received PCM codes. 

The question of digital attenuation PAD mapping has been raised at recent meetings of the' Telephone Industry of 
America (TIA) PCM Modem Ad-hoc Committee. Many participants have been hoping for a single industry standard 
mapping solution. We provide a solution, which will allow PCM modems to function properly throughout the entire tele- 
phone network including with private telephone equipment (PBX's and key systems). 

TIA 464A 

TIA 464A is the industry standard plan for private telephone equipment. Among other recommendations, this spec- 
ification provides the recommended losses between T1 lines and other devices. If a PBX is designed according to the 
loss-plan specification as outlined in sections 4.8.4 and 4.8.5 of 464A, then it should not be necessary to attenuate the 
PCM on incoming T1 lines for on-premise (ONS) connections, since the recommended 3-dB loss is typically imple- 
mented in the analog circuitry of the ONS codec. 

However, Inter -Tel's experience with customers using T1 lines suggests that the recommended 3-dB of insertion 
loss is inadequate and that T1 attenuation pads are indeed necessary. While the standard -3 and -6 dB pads are usually 
adequate, Inter-Tel has provided a wider-range of insertion loss as options for PBX customers. Inter-Tel has imple- 
mented digital gain control on its T1 line cards using two different methods: 

ROM-based mapping 
• , DSP-based algorithm 

ROM Bgsed Mapping 

Inter-Tel first implemented attention pads via an EPROM circuit, which used the incoming PCM code as 8 bits of an 
address and additional address bits selecting the attenuation in 1 dB increments. The output of the EPROM was the 
attenuated digital PCM code. In this implementation, the attenuation range covered from 0 dB to -12 dB in 1 dB steps. 

Experience with customers shows that 0 to -6 dB seems to be the most-often used range of T1 pad values. Taking 
into account the fixed 3-dB of loss included with our ONS card, the net insertion loss for to Digital CO connection is then 
-3dB to -9dB, which mimics well the typical net losses experienced by customers on ONS to Analog CO connections. 

All attenuation tables in the EPROM preserved the LSB of the PCM to preserve the state of the Robbed Bit Sign- 
aling (RBS), if present. However, the preservation of the LSB is unnecessary if no RBS information passes through the 
ROM-based gain control circuitry. 

The ROM-based look-up tables were generated by a custom C program using the following algorithm: 

1. Expand incoming 8-bit,u-Law PCM to corresponding 14-bit, signed integer linear value, x 

2. Calculate output, y, based on equation: y=INT [x*10 (gain/20) )+0.500] 

3. Compress y back into its corresponding 8-bit u-Law value using G.71 1 decision values. 

In addition to performing u-Law to u-Law digital gain control (DGC) through look-up tables, the EPROM also per- 
forms DGC through look-up tables for the following compression schemes: 

u-Law-to-A-law 
A-law-to-u-Law 
A-law-to-A-law 
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DSP Based Mapping 

With the advent of DSP's, for their greater flexibility, the latest implementation uses the Analog Devices 21 xx fixed 
point family of processors. The PCM codes are converted to/from linear values by the ADSP-21xx internal companing 

5 hardware. Inter-Tel's algorithm normalizes the linear value before applying the attenuation multiplier. The attenuation is 
selectable as a linear 1.15 fractional multiplier. The user interface selects the attenuation in fixed, 1 -dB steps over a sim- 
ilar range as the previous implementation. The user configured dB attenuation is converted to the linear fractional mul- 
tiplier by an algorithm. Moreover, the instructions in the Analog Devices DSP can use unbiased rounding, truncation 
and in newer family members, biased rounding. Inter-Tel's implementation uses the unbiased rounding option (RND) 

w during the MAC instruction. This flexibility in rounding and truncation obviously can have a substantial effect on the 
actual attenuation mapping, with the additional caveat of having potentially different mappings for positive and negative 
PCM codes. 

Detailed Description of the Preferred Embodiment 

15 

Detection of the Feature Set 

After a receiver's equalizer has been trained, a known transmit sequence of PCM codes can be sent and received 
using either all or a useful subset of PCM codes. (This sequence may need to be sent multiple times to cover all 6 or 
20 12 bit positions in a robbed bit signaling system.) The receiver can save the linear values for each received PCM code. 
From the known step size between adjacent PCM codes, the receiver can determine those PCM samples that it 
received which are indistinguishable. The indistinguishable samples arise from two transmit PCM codes which when 
scaled by the same attenuation and quantized according to a thresholding procedure create identical PCM codes. 
An example taken from the attached attenuation table (pages A1-A2) which shows the attenuation for each uLaw 
25 PCM code when it is transformed by shifting from 1 to 26 codes numerically. 

Example: 



30 



Transmitted 


Received shifted by 1 PCM code 


uLaw 


Linear 


uLaw 


Threshold 


Attenuation 


dB 


129 


7775 


129 


7647 


0.983537 


-0.14419 


128 


8031 




7903 


0.984062 


-0.13955 



Processing of Indistinguishable PCM Codes 

40 The processing of the feature set relies on knowledge of the coding law (uLaw or Alaw) used. In practice, it is nec- 
essary to consider the preciseness of the numerics used in the attenuation process. Suggestions for accommodating 
implementation deviations is discussed in a later section. The attached C language source code provided as part of this 
application (pages B1- B21) implements a demonstration of the identification procedure indistinguishably codes. The 
basic algorithm is the following: 

45 

1 . For each indistinguishable PCM code received, the original pair of transmitted PCM codes can be determined 
since the PCM codes are sent in a known training sequence. A set of indistinguishable transmitted PCM code pairs 
are determined from this training sequence. 

2. For the first indistinguishable PCM code pair, determine the minimum attenuation required for the larger trans- 
so mitted PCM code to be attenuated to a number of successive lower PCM codes. (A set of 30 successive lower PCM 

codes would be sufficient to cover an attenuation range of approximately 11 dB.) Then determine the maximum 
attenuation permitted for the smaller transmitted PCM code to be attenuated to the same successive lower PCM 
codes. (Minimum attenuation is numerically closer to one, i.e., a lesser attenuation, while a maximum attenuation 
is a greater attenuation.) If the maximum attenuation is a smaller attenuation than the minimum attenuation, discard 
55 the attenuation range as it represents a missing code rather than an indistinguishable code. The minimum and 
maximum attenuations for each successive lower PCM code forms an initial set of candidate attenuation ranges. 

3. For each successive indistinguishable transmit PCM code pair, another set of possible attenuation ranges is. 
determined using the procedure as described in Step 2. Each element of the candidate attenuation range set cre- 
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ated in Step 2 is searched for either full or partial overlap with any element of the possible attenuation range set. 
(In case of partial overlap, the attenuation range in the candidate set can be reduced for refined attenuation accu- 
racy.) In case of no overlap, the candidate attenuation range should be discarded. The discard occurs because two 
different indistinguishable codes must arise from the same linear transformation. 
5 4. Step 3 is repeated for successive pairs of indistinguishable transmit PCM code pairs until either a) the set of can- 

didate attenuation ranges is reduced to a single range or b) all the pairs of indistinguishable transmit PCM code 
pairs have been processed. 

In the first case, the actual attenuation is bounded by the single remaining attenuation range. We can say it is the 
10 detected attenuation can be the average of the high and low attenuations. Alternatively, if the history of each overlap- 
ping attenuation range is maintained (or regenerated) a probabilistic approach of determining a weighted median value 
for determining the detected attenuation can be implemented. 

For the second case where the attenuation appears to remains non-unique, the set of candidate attenuation ranges 
can be used to determine multiple sets of indistinguishable PCM code pairs can be determined for each candidate 
is attenuation. The generated sets of indistinguishable PCM code pairs can then be matched with the original received 
set The generated set using the correct attenuation should match nearly identically if not identically. The sets corre- 
sponding to incorrect potential attenuations would generate significant number of additional pairs of indistinguishable 
PCM code pairs or have missing indistinguishable PCM code pairs. Thus the correct detected attenuation can still be 
determined using this procedure. (Multiple candidate attenuations appear to arise from attenuations of 6 dB or greater 
20 combined with the decreasing step sizes used in with uLaw or Alaw companding.) 

Processing of Missing PCM Codes 

The procedure described above for using indistinguishable PCM codes to determine the attenuation can be 
25 adapted to operate with missing received PCM codes. The transmitted PCM codes corresponding each of the received 
PCM codes around the missing PCM code can be determined. Follow the same algorithm except that the minimum and 
maximum attenuations are exchanged in the algorithms. In other words, the larger transmitted PCM code produces the 
maximum attenuation and the smaller transmitted PCM code produces the minimum attenuation. With this adaptation, 
the algorithm described remains essentially the same. 

30 

Example: 



35 



Transmitted 


Recieved shifted by I PCM code 


uLaw 


Linear 


uLaw 


Threshold 


Attenuation 


dB 


144 


3999 


144 


3935 


0.983996 


-0.14013 


143 


4191 




4063 


0.969458 


-0.26942 



Other Variations 

For robustness with real signals corrupted by noise, it may be desirable to allow a certain number of pairs of indis- 
45 tinguishable PCM codes not to match a candidate range. An alternative to simply dropping a candidate attenuation 
range would be to use a probabilistic greatest likelihood model for selection of the detected attenuation. 

Furthermore, given that implementers of digital attenuation pads may understand and implement the mapping 
process slightly differently (i.e., use of G.71 1 threshold values or average between PCM codes) and allow various inac- 
curacies in the numerical operations (numeric representation, rounding, truncation, etc.), each candidate attenuation 
so range may need to be broadened to make overlap again more likely. The broadening can be by a fixed value or a rela- 
tive amount about the attenuation computed. The previous suggestion of discarding a certain number of pairs of indis- 
tinguishable PCM codes also addresses this problem as very few pairs are likely to be affected by these numerical 
inaccuracies. 

As for digital attenuation pads that pass through the robbed bit signaling (RBS) information, the algorithms pre- 
ss sented can also be applied except that the indistinguishable transmit PCM code pairs will be separated by one code 
and a second indistinguishable transmit PCM code pair will adjacent. The candidate attenuation ranges become bigger, 
but still can be matched by the procedures presented in this invention. 

The procedures described in this invention suggest a robust manner to identify the digital channel attenuation 
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though the use of received PCM codes corresponding to a known transmit PCM code sequence. The use of such an 
algorithm would be necessary when the FCC permits the channel power to be increased in compensation for the sys- 
tematic digital attenuations. 

We further suggest that knowledge of the exact received PCM code set may be preferred by both the analog and 
s digital PCM modems. Echo cancellation in the digital PCM modem may operate with incrementally better performance 
with the knowledge of actual PCM codes presented to the Codec. The receiver's desired spectral shaping can be more 
precisely honored by having the transmitter use the knowledge of the actual PCM codes presented to the Codec in the 
shaping function implementation, then reversely mapping the PCM code after shaping to the transmit PCM code to be 
sent. 

10 The exclusive use of pre-identif ied attenuation table mappings would make the PCM modem technology incapable 
of operating with existing private telephone equipment. It is recommended that the attenuation pad mappings be dis- 
covered for each connection. 

Claims 

15 

1 . A method of determining digital channel attenuation, comprising: 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said attenuation within 
said digital channel; 

20 quantizing said received known training sequence of PCM codes according to a predetermined thresholding 

procedure; 

identifying identical PCM codes created as a result of said thresholding procedure; and, 

determining said attenuation of the digital channel based upon said identification of identical PCM codes. 

25 2. A method of determining digital channel attenuation, comprising: 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said attenuation within 
said digital channel; 

quantizing said received known training sequence of PCM codes according to a predetermined thresholding 
30 procedure; 

identifying PCM codes omitted as a result of said thresholding procedure; and, 

determining said attenuation of the digital channel based upon said identification of omitted PCM codes. 
3. A method of determining a digital channel PCM code transformation comprising: 

35 

receiving a known training sequence of PCM codes, which PCM codes are subjected to said PCM code trans- 
formation within said digital channel; 

quantizing said received known training sequence of PCM codes according to a predetermined thresholding 
procedure; and 

40 determining the transformation of transmitted codes to those received. 

. 4. A method of improved echo cancellation in a communication network having an analog and a digital modem.com- 
prising: 

45 saving codes transmitted from said digital modem to said analog modem for echo cancellation; 

transforming, by a mapping table, codes transmitted from said digital modem to codes received by said analog 
modem; and 

using said received codes as a reference signal for cancellation of echo. 

so 5. A method of improved spectral shaping using a transmit shaping transfer function in a communications network 
having an analog and a digital modem, comprising: 

transforming, by a mapping table, codes transmitted from said digital modem to codes received by said analog 
modem; 

55 using said received codes for transformation to their linear value equivalent representations; and, 

applying said linear value representations to said transmit shaping transfer function. 
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(C) 1997 VoCAL Technologies Ltd 

ALL RIGHTS RESERVED. PROPRIETARY AND CONFIDENTIAL 

VoCAL Technologies Ltd 
3032 Scott Blvd 
Santa Clara. CA 95054 

Product C 

Module. PCM 

This file contaii^s the PCM mapping functions. 

Revision Number- $Revision$ 

Revision Status $State$ 

Last Modified $Date$ 

Identification $ld$ 

Revision History; $Loq$ 
Revision 1 0 1997/03/0) 00:00.00 VD 
Initial release of software 



^include "standard h 
^include "pern h" 
#include -'stdio.h> 
^include <math h> 



r 



struct lavv_array_s { 
sint 1 5 code: 
sint 1 5 linear . 
sint 1 5 other; 
sint 1 5 deselect. 

} 

struct law_s ( 

float tx_poiver_db. 
float rx power _db. 
sinl 1 5 COi int. 
smt 1 13 dnnn 

smt 15 indistinguishable, 
sint !5x, 
sint3 1 bit_iate. 

struct law array s value[T28]. 
} u (aw u (aw rx'~ 

tfdefmeLAW SELECTED 0 
#define LAW DESELECTED MiN 0x0001 
^define LAW DESELECTED MAX 0x0002 
^define LAW DESELEC I ED~DMIN 0x0004 
^define LAWjDESEt EC LED" AVOID 0x0008 
^define L A\AM)E SELECT ED_DUPLICAfE 0x00 10 
#define LAW DESELECTED POWER 0x0020 



^define MAP J"-"RAMF_SiZE 

uint48 frame_sloMnodulus|MAP_FRAME_SIZE]; 
uint48 total__symbols_per_frame 

/ * 

^define USE_4D_TRELUS 

/ * * **** 



void 

u_law_init (struct law_s *law) 
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sir it 15 idx pcm. 

law->tx__power_db = 0 0. 
law--*rx_power_db =0 0: 
law->count = 0. 
Iaw->indistinquishable = 0. 
Iaw->dmin = 0. 
law >bit_rate = 0L. 

for (idx = 0. idx < 128. ♦■ Mrix) { 
pcm - 255 - idx. 



law- > value 
I aw- > value 
law- > value 
taw- > value 



idx 
idx 
idx 
idx 



code = pcm: 

linear = u law_pcm_ decode (pcm), 
other = u~law _pcm threshold (pcm): 
deselect = LAW SELECTED: 



V 



#define SELECT LARGER 



void 



/* appears to work better with standard 3 arid 6 dB attenuation */ 



u_law_attenuate (struct law_s 'law. float attn_db) 

sint 1 5 idx. pcm. 
float attn. 

attn = dbjojloat (-attn_db): 

law->tx_power_db = 00 
law->rx_power_db = 0 0. 
law-count = 0. 
Iaw->indistinauishahle = 0 
law->dmin = t). 
Iaw->bit_rate = 0L. 

for fidx = 0. idx < 128. *■ Hdx) { 
pcm = 255 - idx. 



other = (sint 15) ((float) (u law_pcm_decode (pern)) * attn). 
code = u_law_pcm_encocfe Oaw->vatue[idx]. other), 
linear = u law pcm decode (law->value[idx].code). 
deselect =" LAW SELECTED: 



law->value[idx-l] deselect = LAW_DESELECTED_DUPLICATE: 



law->value[idx 
law->value idx 
law > value idx 
law->value[idx 
// if (idx ■- 0) { 

// if (law->value[idxl code == law->vatue[idx- 1] code) { 

//#ifdef SELECT_LARGER 
// 

//#else 

// law->valiie[idx] deselect = LAW DESELECTED OUPLICATE, 

//#endif 

// } 
" > 



void 

u law range (struct law s *law. S(nl 15 min. sint 15 maxt 



{ 



sint I 5 idx. 

for (idx = 0. idx < 128. * Mdx) { 

if f !aw->value[idx] linear < ruin) { 

law->value[idx). deselect |= LAW_DESELECTED_MIN. 



> 

if (iaw->value(idx) linear > max) { 

law->value[idx]. deselect |= LAW_DESELECTED_MAX; 



void 

u_law__dmin (struct law_s law, sint 15 dmin) 
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sint 1 5 idx. pern; 
sint15 current; 

current = law->vafue[1 25]. linear; 

law->\/alue[ 1271. deselect 1= LAW DESELECTED AVOID; 
law->value| 126]. deselect |=LA\AT[DE SELECT ED "AVOID; 

for (idx = 124; idx >= 0; -idx) { 

if (law->valuelidx|. linear > (current - dmin)) { 

law->value|idxj.desetect |= LAWDESELECTEDDMIN; 

else { 

current = law->value[idx]. linear; 

void 

ujawdiscardjndinstinguishable (struct law s *law) 
sint15 idx; 

for (idx = 1; idx < 125; ++idx) { 

if (law->value[idx] code == law->value[idx-11.code) { 
//#ifdef SELECT LARGER 

// ™ law->valuelidx-1] deselect |= LAW DESELECTED DUPLICATE; 

//#else - ~ 

// law->value(idx].deselect |= LAW DESELECTED DUPLICATE; 

//#endif 

if (law->valueridx-1J.deselect&(LAW DESELECTED DMIN | 
LAW DESELECTED AVOID | LAW DESELECTED POWER)) { 

law->vatue[idx-1] deselect |= LAW DESELECTED_DUPLICATE; 

else { 

law->value[idx) deselect |= LAW_DESELECTED_DUPLICATE; 

, - 1 ' 



************************* 



void 

u_law_attenuate_dmin (struct law_s Mxjaw, struct law_s *rxjaw) 
sint15 idx; 

rx law->value[1 271. deselect 1= LAW DESELECTED AVOID; 

3] i ■ 



rx~law-> value] 126] deselect |= LAW^DESELECTED^AVOID; 

for (idx = 124; idx >= 0; -idx) { 

ifftx law->value]idx] deselect & LAW DESELECTED DMIf 
rx_law->value[idx] deselect |=~LAW DESELECTED 



t : 

void 

ujaw exclude (struct law s *law) 

sint15 idx, pem; 

for (idx = 127; idx >= 0; --idx) { 
pem = 255 - idx; 

if ((pem == 165) [| (pem == 169)) { 

^ law->value[idx) deselect |= LAW_DESELECTED_AVOID; 

r *************************** ****.************* ********.*** v 

float 

u_law_power (struct law_s *law, float max_db) 
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sint 1 5 idx. count, dmin. prev; 

float factor, power, max_float; 

float tx_power. tx_sum, rx_power, rx_sum; 

factor = pow(10.0. ((3.17 + 3.01244) / 20.0)) / 8159.0; • 

count = 0; 
tx_sum = 0.0; 
rx_sum - 0.0; 
dmin = 8159: 
prev = -8159; 

max float = db_to power (max_db); 
for (fax = 0; idx < f28; ++idx) { 

if (law->value[idxj. deselect == 0) { 

rx_power = law->value[idx|. linear * factor; 

rx_power = rx_power * rx_power; 

rx_power = rx_power + rx_sum; 

tx__power = ujaw.value[idx]. linear * factor; 
tx_power = tx_power * tx_power; 
tx_power = tx_power + tx_sum; 

power = tx_power / (float) (count + 1); 
if (power <= max__ftoat) { 

rx_sum = rx_power; 

tx_sum - tx_power; 

count+ + ; 

if ((law->value[idx]. linear - prev) < dmin) { 

dmin = law->value(idx]. linear - prev; 

prev = law->value[idx) linear; 

else { 

law->value[idxj. deselect = LAW_DESELECTED_POWER; 

//printf ( ,,0 / 0 f %f '%d %d\n'\ power. tx_sum. count, law->value[idx). deselect); 

!x_power = tx_sum / (float) count; 
tx_power = power_to_db (tx_power); 

rx^power = rx_sum / (float) count; 
rx_power - power_to_db (rx_power); 

law->tx_power_db = tx_power; 
law->rx_power_db = rx_power; 
law->count = count; 
law->dmin = dmin; 



} 

/* 



return power; 



******** 



void 

uJaw_counMndinstinguishable (struct law_s *law) 
sint 15 idx, count; 
count = 0. 

for (idx = 1 ; idx < 1 28; ++idx) { 

- if (law->value[idx].deselect == LAW_DESELECTED_DUPLICATE) { 
count+ + ; 

law->indistinguishable = count; 



/* 



********** 



a********************************** */ 



sint31 

u__law_bit_rate (struct law_s *law) 
float bit rate float; 



bit_rate float = (log 10 ((float) law->count * 2.0) / log 10(2.0)) * 8000.0; 
law->biCrate = (sint31) bit_rate_float; 
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} 

r 



return law->bit rate: 



7 



//#define MAX SHIFT 20 
#define MAX SHIFT 30 
//#define MAX_SHIFT 40 

float 



/* covers to -7.08 dB 7 
r covers to -1 1.2 dB 7 

/* covers to -14 65 dB 7 



ulaw_detect_attenuation (struct law_s Maw) 



{ 

stntIS idx. idx I. idx2. idx3. trim; 
sintl 5 count; 

float attn_low[MAX SHIFT], altn highlMAX SHIFT] 
float a lowiMAX STHIFT). a highTMAX SHIFT) 
float afln. cfiff; "~ ~ ~ ' 

count = 0; 
idx3 = 0; 

for (idx = 128; idx >= MAX SHIFT; -idx) { 

if (law->value[idx) deselect & LAW DESELECTED DUPLICATE) ( 
idx3<- + ; ~ - 71 

if (count == 0) { 

ttfdef SELECT_LARGER ^ (C ° Unt = ° C ° Unt < MAX - SH,FT : -count) { 

ujaw.value[idx + 1| linear; attn_high[count] = (float) ujaw.value[idx-count].other / (float) 

u law. value[idx]. linear: 
#else 

u law value[idx| linear; 

u law. valuefidx- 1 1 linear 
#endif 
// 

attn high[count]. attn Jowfcounty 

count = MAX SHIFT; 

} 



attnJow[count) = (float) u_law value[idx-count-1]. other / (float) 

attn_high[count| = (float) uJaw.value[idx-coiint- 1 ]. other / (float) 
attnjow(count) = (float) u Jaw value[idx-count-2|. other / (float) 

printf ("original %d %f %f\n" ( u_law.vaIuefidxJ.code. 



else { 

#ifdef SELECT_LARGER 

ujaw.value[idx + 1 J. linear; 

ujaw value[idx|. linear; 
#else 

u_law.value[idxj. linear; 

u law. value[idx-1 1. linear 
#endif 

// 

a_high[idx1J. aJow[idxl|); 



a_high[idx2))) { 
attn_high[idx1], attnlow[idx1 1. a_high(idx2]) l 



aJow[idx2|)) { 

attn highfidx 1 1, attnfow[idx1], aJow[idx2|); 



for (idx I = 0; idx 1 < MAX_SHIFT; ++idx1) { 

a_high[idx1) = (float) ujaw. value[idx-idx 1 ] other / (float) 
ajowpdxlj = (float) ujaw value[fdx-idx 1 -1 ]. other / (float) 

a Jiigh[idx 1 ] = (float) ujaw.value[idx-idx I - 1 J. other / (float) 
a_low[idx1J = (float) ujaw.value[idx-idx1 -2] other / (float) 

printf ("generating %d %f %f\n". ujaw. value(idx). code, 

} 

for (idx I = 0, idx 1 < count:* + +idx1) { 
trim = 1 ; 

for (idx2 = 0; idx2 < MAX SHIFT- + Hdx2) { 

if (<attnjiigh[idxfj> = a_high(idx2)) && (attnjowjidx 1] 



> 



printf ("reduce high %f %f %f\n", 
a«n__hg[h[idx1] = a Jiigh[idx2); 



if ((attnjiigh[idx1) >= aJow(idx2]) && (attn low[idx 1 1 < = 
printf ("increase low %f %f %ftn". 



} 



attnJow[idxl] = aJow(idx2J: 
trim = 0; 
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a Jow(idx2))) { 
attn_iov;[idx 1 j). 



EP0 871 303 A2 

if ((3ltn_:ug'.i[idx1] <= o_high!«<ix2!) (rttnjo- v|idx I ] 
rrintf ("hold ^; %fVT. attnjitgnfidx 1 1. 
trim = 0; 



if f'attn_high[idx i J < attn_lo\v[idx I )) { 
tnm = I . 

If Mum) { 

printf ("trim %f %An" attnjirghjidx 1 1. attn low(idxlj): 
for Hdx2 = idxl: idx2 < {count - r ' ~* " 

attnjiigh[idx2j = attn hie 

attn Jo w[ idx 2] = attnTovvji 

--count: 



micix I j. aun i< 
1), +»idx2){" 
iiqhlidx2* (]. 
>w|idx2*- I); 



} 



} 



} 



if (count <= 0) { 

pnntf ('count break\n"l 
break. 

}\ (tdx3 >= 10) break 



IN 

IN 

rangeArV') 



tf (count ^ \) { 

attn = (attn hiqh[0| ♦ attnjow[0]) / 2.0. 

diff = atti>_mgh(Oj - atlnjovv(0j 

for "(idx 1 = 1 .^idx f < count. + -+tdx 1) { 

if ((attn Jiighfidx 1] - attn low[idx1|) < diff) { 

attn = (attn hjghlidxfj f altnjowlidv I ]) / 2 0. 
diff - aUn_fi"igh[tox 1) - attn_low[icfx 1 ). 

pnntf (' Muplitple potential attenuations detected, resolve by generating and malchingMV 
"duplicated codes to those observed Currently choosing largest 



else { 



} 

r 



attn = (altn_high[Oj + attnjo\v[0]) / 2 0 • 
attn = tog 10 (attn) * 20 0- " 

printf ("count "'id attempts %d range %t 'enn'* count. idx3. attn_higb(0). attn !ow|0]j. 
printf ( "attenuation detected %f>n". attn). 
return attn 



sintlS 

u law niap symbol set {stiuct taw_s *tx_taw struct law s *rx_law) 

sint 1 5 fdx count limit 
uint 1 6 mask 

total symbols _pef_hame Isv/ = 1. 
lotal_symbols per frame mid = 0. 
total symbols p^r frame msw = 0. 

// prinlff "%04x °.'o0'1x %0'Lx\n". totai_symbols perftame msw 

// total_symbols_per_(tame mid. toTai_symbols_per_frame Isw); 

for (idx = 0. idx < MAP_ FRAME SIZE: + + idx) { 



frame_$lot_modulus 
frame _s lot modulus 
frame slot" niodulus 



idx 
idx 
idx 



isw - total _symbols_per_frame Isw; 
mid = totaf symbols_per_fTame.mid: 
msw = lota"T_symbcls_per_frame msw: 



#ifdef USG_4D TRELLIS 

"limit = tx !aw->count 
if (idx & T) { 

limit = limit >> 1 . 

} 



B6 



13 



BNSDOCID: <EP 0871 303 A2J_> 



#else 



#endif 

II 
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for (count = 1 , count < limit: + + counU { 

umt'l8_add (&tota!_symboIs per_frame. &frame_slot_modu!us[idx]. 

&to(al_symbols_pet_frame). 

for (count = 1, count < tx_law->count. count) { 

uin(48_ acid <&total_symbols per_frame. &Uat ne_slot_modulus(idx). 
^ &tofal_symbols_per_frame); 

pnntff'%04x %CMx %04x\n" total_svmbo1s per_frame msw. 
j i ^ lolal_symbols__pe»_frame mid toral_symbols_per_frame Isw). 

count = 47 

mask = 0x8000. 

for (idx = 0. idx < 16 *- Mdx) / 

if (mask & totnl_symbols per frame msw) { 
return count; 

/ 

mask = mask >> ) . 
-count; 

} 

mask = 0x8'J'J0: 

for (idx = 0. idx < 16; + Hdx) I 

if (mask & total_symbots_per frame mid) { 
return count. 

mask - mask >> I*. 
-count. 

} 

mask = 0x8000. 

for (idx = 0, idx < 16. + +idx) { 

if (mask S total_symools per frame Isvv) { 
return count. 

} 

mask - mask >> 1 : 
-count; 

return count 

) 

/ / 

void 

crinldeselect (suit 15 deselect) 

if (deselect & LAW DESELECTED_MIN) { 
^ printf ("(Below Mm) "). 

if (deselect & LAW_OESELECTED__MAX) { 
prinlf ("(Above Max) "). 

) 

if (deselect & LAW_DESELECTED_AVOID) { 
printf {"< Avoid Omin) "); 

} 

if (deselect X LAW DFSELFC7 ED__DMIN.) { 
printf ("(Fiefween Dmin) "); 



if (deselect & LAW_DESELECTFD DUPLICATE) { 
pnntf f (Duplicated Code) "T. 



if (deselect & LAW_OESELEC F EU POWER) { 
printf ("(Power Limit) "), 

} 

prinff("*.n"). 
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//^define DO J : EW 
/ *" 



//^define UO_PRINT 0 
//^define DO PR1NT~1 
//^define DO~PRlNT"2 
//^define DO~PKINT"3 
//^define DO 'PRINT 4 
Wdefine DO~PRINT~5 



#ifdef DO l"EW 
fldefine MAX AH N ) 
#define MAXDMIN 4 

#define START DMIN 48- 
tfdefme SI FP DMIN 4 

//define STAR T AH N 110 
^define STEP ATTN I o 



#clse r DO FEW V 
#dofine MAX ATTN 14 
itticUnc MAX" DMIN 36 

^define START DM/N 4 
^define S T FP_DMIN 4 

tfdefine START ATTN 0 0 
fldefine STEP ATT N 1 0 
tfcndif r DO PEW 7 



void 

main (void) 

{ . 

stnt15 idx. count: 
float attn 

sint 1 5 bits per_frame: 

sint31 bit rate. overal!_bit rate: 

float input attn 

S" 4 1 ft inpt it tx_dmm input ix drmn 
smt 1 5 i. j. results idx. 

struct resiilts_s { 

float attn. 

sint t 5 drnin. 

sint31 bit_rate_simple; 

sint3 1 bit rate exclude. 

sint31 bit" rate /distmguis.habl9_tx_dmin. 

c .int31 bit rate distinquishable_rx_dmm 
} results (MAX ATTN r MAX JDMIN] 

results idx - 0. 

input._attn = STAR r ATTN. 
// • inpuMx_dnnn = 12CT 
// »npi/t_rx_dniin = 120; 

for (i = 0: i < MAX ATTN, input attn += STEP_ATTN. i + + l { 
inpuMx'jTmin = S 1 ART DMIN: 
input Vx dmin = START~DMIN, 

for (j = 0. j < MAX_ DMIN. inpuMx_dmin += STEP_DMIN. input_rx_dmin += STEP_DMIN. j + + ) { 

p/intf ("\n"). 
u_law_init (&u_.1aw); 

// ujaw_range (&u law. 30. 4000). 

ujawjdmm (4u_Faw. input tx dmin) 



// u_law_pcwer (Sujaw. -12.D C-Q 47981V): 

u law_power (&u_law. - 1 2.0): 
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■ ujaw__bit_rate (&u_law). 

pnntf C°/c6 21 °o3d %3d ". input_attn. input Jx_/Jmin. input_rx_dmin). 

fl printf ("tx power ° 0 3 2f rx power °c3.2f count °:>d dmin %d dup \cd bit rate %ld\u" 

printf ("tx pwr 3 2f ix pwr %3 2f cut %2d dmin %3d dup %2d rate %ld\n". 

ujaw tx_ power db. ujaw ;x_power db. ujaw count, ujaw dmin. 
u law. indistinguishable, ujaw bit rale); . 

#ifdef COMPRINT 0 

for lidx = 0. idx < 128. ^id*H 

printf ("% 3d %5d 'vobd %5d %5d ujaw. value[idx) code, u law. value[idxl. linear . 

u law valuejidx) other u_law pern encodefu law value(idx) linear ). 
u law pem encodem Jaw. valuefidxj other - fj). 
puntdeselechu law vaiuefidx] deselect). 

#endif * 



u law_atteiuiate (&ujaw rx. input attn); 
11 ulawrange (8-u Jaw_rx~22. 40007. 

// tflaw clrntn (&u law_rx. input_rx dmin). 

ujaw attenuate" dmin (&u law. Sujaw rxi /* transfer from transmitter "/ 
// ujaw_exclude \SuJaw jxf. 

U u_law_discard indistinguishable (&u law_rx). 

u law power (~&u law rx. -12 0); 

Li law discard incIinstlhgLJistiabJe (&ujaw_rx). 

u_ lavv_count indinstinquishable (&u Jaw_rx). 

ujaw_bit rafe (&u lavv rx); 



printf 2f %3d c /n3d ". inputattn. input tx /imtn. input_rx dmin). 

pnntf ("tx power ' ; .\3 2f rx power %3 2f count %d dmin %d dup %d bit rate %ld\n" 
printf ("tx pvvr %3 2f rx pwr %3 2f cnt %2d dmin %3d dup %2d rate %ld\n'\ 

u_ law_rx tx_power_db. u taw rx rx_power_db. ujaw_jx count. 

ujaw rx dmin, u law_rx.TndisTingiiis liable. u_law rx.bit rate). 

#ifdef DO_PRlNT 1 

for (idx = 0. idx < 128. + + tdx)J 

printf ("%3d %5d ":obd %5d ". u_ law value[idx] code, ujaw rx value[idx] other, 
u law rx value[idx]. linear . ujaw^rx.valuelidxj! code); 
^ print_deselecl(ujaw_rx valuefidx] deselect). 



#endif 



attn = u law_detect attenuation (&u iaw_rx). 
printf ("attenuation detected %f\n". altnV 

resnitsf results _idxj nt'n -input attn 
resultsfresults Jdx] dmin = inpuT_tx_dmin: 
if (ujaw_rx dmin == 0) { 

(csults[results_idxj bit_rate_simple = 0. 

else { 

. results[results jdx) bit_rate_simple = ujaw_rx,bit_rate: 



ujaw attenuate t&ujaw rx. input attn): 
// u_ lawjrange (&u law_rx.~22. 40007. 

// ujaw_dmin (&Li_raw_rx, input_rx dmin). 

u_law_attenuate dmin (Su law, £*ujaw_rx). /' transfer from transmitter */ 

u law exclude (5u faw_rx^ 
// u Jaw__discard indistinguishable (&ujaw_rx) 

ujaw" power f&u law^rx. -12.0); 

u Jaw_dtscauJ_inOinstmyuishabie (&u Jaw_r v ) 

u_law__count indistinguishable (&ujaw_rx). 

u Jaw_bit_rafe (&ujaw_rx). 

printf ("Vo6.2f %3d %3d ". iput_attn. ipuMx_dmin. input_rx_dmin): 

// printf ("tx power %3.2f rx power %3.2f count %d dmin %d dup °/od bit rate %ld\n". 

printf ("tx pwr %3 2f rx pwr %3.2f cnt %2d dmin %3d dup %2d rate %ld\n". 

u_law_rx tx^power db. u jaw rx.rx nower_db. u_law rx count. 
ujaw_rx dmin. ulawrx indisTmguisnable. ujawrx.bitrate): 

#ifdef DO PRINT 2 
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for (idx = 0. idx < 128. *-Hdx){ 

pnnti ("%3d °o5d %5d %5d ". u law valuefidx], code, ujaw rx.vnlue(idx). other. 

u law_rx.value[idx]. linear. u_law_rx value(idxlxode); 
print_deselect(u_faw_rx.value[idx] deselect): 



tfendif 



atln - u lav/ detect attenuation (&u law_rx) 
printf ("attenuation detected %f\n". altn); 



if Cu_lavv_rx dmin == 0) { 

results[results_idx| bit_rate_exclude = 0. 

} 

else { 

results[results_idx] bit__rate_exc(ude = u Jaw_rx.bit_rate: 



ti low_attenuate (&u_law rx. input attn): 
u"law_rangc {^u law rx. "22. 40007. 
if law_dmin (&u_Taw_rx input_rx dmin); 
if law attenuate_dmin (&u_law. <STu law_rx). 
ujaw discard indistinguishable (Su law_rx>: 
1 1_ law po we r ( S 1 1 I a vv_ r x . - 1 2 . 0) : 
ujawjpcwer (&ujaw_rx. (-12 0 - attn)). 
// uJaw_count indistinguishable (&ujaw_rx). 

u law bit_rafe (&u lawrx) 



7 



printf 2f %3d %3d ". fnpnt_attn input Jx_dmin. inputrxdmin); 

// printf f'tx power n ;-3 2f rx power %3 2f count %d dmin %d dup ?cd bit rate %!d^n". 

pnnif f tx pv/r r '-3 2f rx pwr %3 2f cnt %2d dmin v«3d dup %2d rate %ld\n". 

u law_rx tx_power_db ujaw rx rx_ powered b. ujaw_rx count 
ujaw_rx dmin. ujaw_rx indisTinguisnable. li_ law_rx.bit_rate); 

#ifdef DO_PRfNr 3 

for (idx = 0 idx < 128. + Hdx) { 

pnnif f" J ?3d %5d />5d %5d " u law valife(idx) code, ujaw rx.value(idx) other. 

u Iaw_rx.value[idx7. linear, ujaw /x.vatue|idxf code); 
pnnt_deselect(u_Taw_rx vatuefidx] deselect), 



ffendif 



r 



attn = ujaw_detect attenuation (&u law_rx). 
printf ("attenuation detected %f\n" attn), 

iesults[results_ tdx) bit_r ate_distinguishable_ tx chuin = ujawrx bit_rale. 



u law_attenuate (&ujaw rx. input attn), 
// ujaw^rancjc (&u lawj-x.~22. 4000X 

u_law_dmin (&u_Taw rx. input_rx dmin). 
// n_law_attenuate_dmi7i (&ujaw. <5u law_rx); 

ujawdiscard indinstinguishable (&ujaw_rx). 

u law_ power ("Sujaw rx. -12 0). 
// u"tnw power (3u_law~_>x. (- 12 0 attn)). 

// uJaw"count indnistinguishabie (&ujawjx) 

1 1 Jaw " hit r a fe ( F. \ i la w r x ) . 

printf 6 2\ c "3d %3d ". iuputattn. tnpulJx_rUntn. input _rx_<imtn) : 

// printf ("|x power %3 2f rx power %3 2f count %d dmin %d dup %d bit rate ^ld\n". 

printf ("tx pwr %3 2f rx pwr %3.2f cnt %2d dmin %3d dup %2d rate %ld\n", 

u_taw_rx tx_power_db. ujaw rx rx_power_db. ujaw_rx count. 
ujaw_rx dmin. u_!aw_ rx indisTinguisnable. u law rx bitjate). 

#ifdef DO_PRINF 4 

for (idx - 0 idx < 128. +<idxW 

pm»tf C%3d %5d %5d %5d ". u law. vaiue(idx). code, ujaw rx. value(idx) other, 
u law_rx.value(idx], linear, u Jaw_rx.value[idx]xode): 
^ phnt_deselecKu_Iaw_/x.value[idxJ deselect); 

attn = u fawdetect attenuation (&u law_rx); 
printf ("altenuation cfetected %ftn" ; altn); 

^endif 
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resiJlts|resul!s_idx + +] bjt_rate_distingipshab!e_rx_dmin = ujaw_rx. bit_rate. 

' / 

bits_per_ffr»me = uJaw_map_symbol_set (&u_law_rx. &u_law_rx): 

bit rate = (stnt3 1 ) (((float) (bits_per_frame) / 6 0 + t 0) 4 8000 0): 

#ifdef USE_4D TRELLIS 

overall~M_ratc = (sint31) ({(float) <bits_per_frame) / 6.0 +10 + 0.5) * 8000 0)' 

#else 

overall._bir_rate = <sint31) (((float) <bits_per frame) /6.0 +10)* 8000 0). 

#endif 

#ifdef DO_PRINT_5 

pnntf ("bits per frame 3 od bit rate %\ti overall bit rate %ld (with trellis)\rV\ bils_per frame, 
bit rote, overall bit rate): 
#endif - _ - 

#ifO 

for (count = 0. count < 1000. count = count + bits_per_frame) { 
gene'ate_bits (bits_per_frame). 
u Jaw_send_ fiame (&u iaw &u law__rx). 
u la\v_receive_frame (5u_law. 5Tu_law_rx). 
cFTeck__bits <bits_per_frame). 

#endif 



tfifndef DO rEW 

for^idx = 0; idx < MAX ATTN * MAX DMIN. idx^ + ) { 

printf("%6 2f "?-Qd %5ld ^SIcHoSId %5ld\r/ resultsfidxl.attn. results[idx].dmm. 
" 11 1 bit_rate_simple. results(idx) bit_rate_excluae. 
bit_rate_distinguistiabte_tx_dmin, 
bit_rate_distinguishab1e_rx_dmin); 



#endif * 



results 
results 
results 



idx 
idx 
idx 
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pcmcnv.c 

(C) 1 997 VoCAL Technologies Ltd 

ALL RIGHTS RESERVED PROPRIETARY AND CONFIDENTIAL. 

VoCAL Technologies Ltd. 
3032 Scott Blvd 
Santa Clara. CA 95054 

Product: C 

Module PCM 

This file contains the PCM conversion functions. 



Revision Number: 
Revision Status 
Last Modified 
Identification: 



$Revision$ 
$State$ 
$Uate$ 
$ld$ 



Revision History: $Log$ 
Revision 1 0 1997/03/01 00:00:00 VD 
Initial release of software 



//include "standard h" 
//include "pcm h" 



//include <stdio h> 
sint 1 5 

u_!aw_pcm_encode (sfrac! 15 input) 
sfractl 5 value, seq; 



value = input; 

if (input < 0) value 



-input. 



value = value + 33; 

if (value > 8159) value = 8159. 

seq = iexp (value); 
//printf("%5d %5d %5d", input, value, seg): 
value = norm (value, seg): 

value = value A 0x4000; 
value = value >> 10; 
//printff %5d", value): 

if (input < Oi { 

value = value | 0x0080. 



seg 
seg 



seg + 9; 
seg « 4; 



value = value i seq; 
lue); 



//printff %04x M . vali 



value = value A OxOOff; 
//printff" %5d\n". value): 

return (sint 15) value: 

} 

#ifdef DO_U_LAW 
x 
x 

X 
X 
X 
X 
X 



/* Take absolute value 7 

/* Limit values */ 

/* Compute exponent */ 

/* Normalize value 7 

r Clear bit V 
/* Position bits 7 

/* Insert sign bit 7 

/* Compute segment 7 
/* Insert segment bits *f 
/* Invert bits 7 



AR = ABS AR; 
AY0 = 132; 
AR = AR + AY0; 
SRI = 32636; 
IF AV AR = PASS SR1. 
SE = EXP AR (HI), 



r Input sample passed in AR (1.15) 7 
/* Take absolute value 7 
/* 33 « 2 7 

/* Limit values 7 

- /* Compute exponent 7 
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y 
x 
x 

X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

tfendif 



f Normalize value 7 



SR = NORM AR <LO): 
AYO = 0x4000; 
AR = SRO XOR AYO. 
SR = LSHITT AR" BY -10 (LO): 
AXO = SF. AR = PASS AYO. 
ir POS AR = PASS 0. 
SR = SR OR LSHIFT AR BY -7 (LO): 
AYO = 7. 

AR = AXO *- AYO. 

SR = SR OR LSHIFT AR BY 4 (LO) 
AYO = 0xFF 
AR = SRO XOR AYO: /* Invert bits */ 

RTS. /* Pass on rate converted value 



r Clear bit 7 
r Position bits 7 

/* Create siqn foil 7 

/* Insert sign bit 7 
r 9 if using 803 1 range 7 
/* Compute segment 7 
/* Insert segment bits 7 



/* 



sfract 15 

u_law_pcm_decorte (umis code) 
sfract t5 value, seg 



value 
value 



code & OxOOff: 
value A OxOOff: 



sea - value & 0x0070: 
vafue = value ' i seg 
seg ~ seg >> 4- 

if ((value - 0x0080) >= 0) { 



Mask bits 7 
Invert bits 7 

Isolate segment bits 7 
Remove segment bits 7 
r Shift to tow bits 7 



else { 



} 



value 
value 
value 
value 



value 
value 
value 



value » Oxff80 
value value + 33 
value << seg: 
33 - value: 



value + value + 33. 
value « seg. 
value - 33; 



r Determine sign 7 



return value. 



Process negative value 7 
/* Add seqment offset 7 
(* Position bits 7 

/* Remove segment offset 7 

/* Process positive value 7 
/* Add segment offset 7 
/* Position bits 7 

r Remove segment offset 7 



sfract 1 5 

u law pem threshold (uintS code) 

{~ ~ " 

sfract 15 value, seg. 

value = code & OxOOff; 
value = value OxOOff. 



r Mask bits 7 
r Invert bits 7 



seg = value & 0x0070. 
vafue = value S OxOOOf 
seg = seg >> 4 
seg ^ ♦ 

value = value * 1 7. 
value = value < -' seq. 
value = value - 33: 

if ((code S 0x80) == 0) { 
value = -value 



> 

return value 



r Isolate segment bits 7 
/' Remove segment bits 7 
r Shift to low bits 7 



r Add segment offset 7 
/' Position bits 7 " 

/* Remove segment offset 7 

/* Determine sign 7 

/* Process negative value 7 



#ifdef DO__A_LAW 
aJavv_pcm__encode 



AR = ABS AR. 

AYO = 511: /' 127 7 ■ " 
AF = AR - AYO: 

IF GT JUMP ajaw_pcm_enc_l 



r Input sample passed in AR (1.15) 7 
/* Take absolute value 7 
/* Check for zero segment 7 
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SR = LSHIFT AR BY -4 (LO). r Downshift bits 7 
AR = 0x4000. 

IF NEG AR = PASS 0. /* Create sign bit 7 
SR = SR OR LSHIFT AR BY -7 (LO); /* Insert sign bit 7 
AYO = 0x55: 

AR = SRO XOR AYO /* Invert bits 7 

RTS. /* Pass on rate converted value 7 

SE = EXP AR (HI). r Compute exponent 7 
AXO = SE. SR "= NORM AR (LO) /* Normalize value 7 

AYO = 0x4000. 

AR = SRO XOR AYO, /* Clear bit 7 

SR = LSHIFT AR BY -10 (LO). r Position bits 7 
AR - PASS AYO. 

IF NEG AR = PASS 0. /* Create siqn bit 7 

SR = SR OR LSHITT AR BY -7 (LO); /* Insert sign bit 7 

AYO =7. 

AR = AXO + AYO; /* Compute segment 7 

If I ! AR = PASS 0: 

SR = SR OR LSHIFT AR BY 4 (LO); /* Insert segment bits 7 

AYO = 0x55. . 

AR = SRO XOR AYO. r Invert bits 7 

RTS. r Pass on rate convened value 7 

idenKIDENt RX B1 PCM DECODE). 

a_law_pcm_ decode Input sample passed in AR 7' 

AYO = 0x55 /* Ma^k bite for inversion V 

AR = AR XOR AYO: r Invert bits 7 

SRI = 0x3, /* Set sign bit 7 

SRO - 0x800: r Set LSB.of interval 7 

SR =*SR OR LSHIFT AR BY -12 (LO). /* Isolate seqment and interval 7 

AYO - 32. 

AXO = AR. AF = PASS AYO. 

AYO = 9; . r Segment bias 7 

AR = SRI - AYO, r Determine shift 7 

IF LT AF = PASS 0, r No extra MSB bit 7 

IF EQ AR = PASS.O. r No less than zero bits 7 

SR = LSHIFT SRO BY - It (LO): r Isolate interval 7 

SE = AR. AR = SRO OR AF: /* Add bit if necessary 7 

SR " I SHIFT AR (LO). /* Position output 7 

SR = LSHIFT SRO BY 3 (LO): 

AY0 = 0xrF80 

AR = SRO. AF = AXO + AYO: /* Check if sign bit set 7 

IF LT AR = - SRO. /* if it is. negate result 7 

RTS /* Pass on rate converted value (I 15) 7 

#cndif r DO A" LAW 7 
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* * / 

math c 

(C) 1997 VoCAL Technologies Ltd. 

ALL RIGHTS RESERVED PROPRIETARY AND CONFIDENTIAL 

VoCAL Technologies Lid 
3032 Scott Blvd 
Santa Clara. CA 95054 

Product C 

Module MATH 

this file contains fractional math support functions 

Revision Number SRevisionS 

Revision Status $State$ 

Last Modified SDateS 

Identification $ld$ 

Revision \ listory. SLoqS 

Revision 1 0 1997/03/0 TOO 00 00 VD 

Initial release of software 



tfincti ide "standard h" 
^include "pcm.h" 
^include <ma!h h> 



sint 1 5 

iexp (sfract 1 5 value) 
{ 

sint 1 5 exp. 

if (value < 0) { 

value = -value. 

} 

exp - } 5 

if (value = = 0) { 
return exp 



exp - 0, 

while ((value & 0x4000) == 0) { 
-exp. 

value = value << 1 . 

} 

return e*p 



sint 1 5 

nor m fsfract I 5 value sint 1 5 exp) 
{ 

return (value « -exo* 

) 



r 



float 

dbjojloat (float db) 

^ float db_float. 

db float = db / 20.0. 

dbjloat = pow (10 0 dbjloat). 

return db float; 
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} 

/' 



float 

db to power (float db) 

' fioat dbjloat: 

dbjloat = db / 10 0. 

dbjloat = pow ( 10 0. dbjloat). 



} 



return db float 



float 

floaMo_rib (float power ) 
* float float_db. 

floaPdb = log 10 (power) • 20 0: 

return float db. 



float 

oower_to_db (float power) 
float powered b. 

powerdb = log 10 (power) * 10 0: 
return power db. 



void 

uint48_add (uint48 *srd. uint48 *sk;2. uinM8 *dst) 

uint 1 6 s 1 1, s 1 in. s 1b. s2l s2m.s2b.dstt.dstrn.dslb. 

s 1 1 = src 1 ->lsw. 
s Ini = src 1 ->mid. 
sib" 3 src 1 ^msw; 

s2l " src2->lsw. 
s2rn = src2 >mid. 
s2b = src2- :> msw: 

asm { 

mov sir 
add ax. s2l, 
mov dstl ay. 

mov av s 1 in. 
adc ax. s2m, 
mov dstm. ax. 

mov ax s lb. 
adc ax. s2b. 
mov dstb ax. 
} " 

dst->!sw = dstl; 
dst-^mid = dstm. 
dst->msw - dstb. 

// printf("%04x %04x %04x - %04x %04x %04x + %04x %04x %04x\n" 
// dstb dstm. dstl. slh. slm. sH. s2h. s2m. s2l);. 

) 



void 
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uint48_sub (uint4S *srcl. uint48 *src2. uinMS *dst) 

uinM6 s II. slm. s1h. s2l. s2m. s2h. dstl. dstm. dsth: 

sll = src1->lsvv. 
slm - src 1 ->mid. 
s Ih = src I ->msw 

S2I = src2->lsw. 
s2m = src2->niid. 
s2h = src2->msw. 

asm { 

mov ax s 1 1. 

sub ax. s2l. 
mov dstl ax. 

fiiov ax s I m; 
sbb ax s2m. 
mov dstm. ax. 

mov ax. s lb 
sbb ax s2h. 
mov dsth ax. 



dst->Isw = dstl: 
dst->mid = dstm 
dst->msvv = dsth 



sfract 15 

frnpy (sfract 1 5 a. sfract 15 b) 

long a lonq, b long, 
short a_shbrt 

ajong - a, 
b_long = b. 

ajong = ajong * bjong: 
a_short - (short) (a Jong >> 15): 

return a short. 

} 
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r 

pcm h 

; fC) I997 VoCAL Technologies Ltd 

ALL RIGHTS RESERVED- PROPRIETARY AND CONFIDENTIAL 

* VoCAL Technologies Ltd 
3032 Scott Bivd 
; Santa Clara. CA 95054 

; Product: C 

Module PCM 

This file defines the PCM utility functions 

Revision Number $Revision$ 

Revision Status SStateS 

Last Modified $Date$ 

Identification StdS- 

Revision History SLog$. 



7 

r * 



#ifndef PCM PCM H 
^define PCM PCM M 



typedef struct ( 

uint 1 6 Isw, 
uint t6 mid. 
uint 16 msw: 

} uint48: 

sint15 iexp (sfract 15 value) 

smt 1 5 norm (sfract 1 5 value, sint 1 5 exp): 

float dbjojloat (float db). 

float dbjo power (float db); 

float float_fo db (float power), 

float power_fb_db (float power) 

void uint48 add (uint43 'srcl. uint48 *src2. uint48 *dst) 
void wint48_suh (uint48 'srct. uint48 *src2, uint48 *dst): 

stnt 15 u law_pcm_enccde (sfract 15 input), 
sfract I5~ii_law_ pcm_decode (uints code): 
sfract 15 u law "pcm threshold (uint8 code) 

r * ...... . • — .. — ...... 

#endif /' PCM PCM H 7 
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/ * * * k ----- v 

r 

\ standard h 

\ (C) 1994. 1995 1996 VoCAL Technologies Ltd. 

; ALL RIG* ITS RESFRVF.D . PROPRIETARY AND CONFIDENTIAL. 

' VoCAL Technologies Ltd. 
3032 .Scott Blvd 
\ Santa Clara CA 95054 

; Product MODEM 101 

; Module SYSTEM 

This file defines the standard system definitions 

Revision Number $Revision$ 

Revision Status SStateS 

Last Modified $Date$ 

; Identification. $ld$ 

Revision History SLogS 

*/ 



tfifndef SYSTEM STANDARD H 
#dcf:nc~ SYSTEM" S TANDARfJ" H 



/* 



#ifdef DO_CX025 
include <GLBL Clonals h> 
//include <m68EC040 h> 
#end:f r DO_CX025 7 

//ifdef DO_W1N 

/'-//incfude < null h>-7 
#define STDC FALSE 
/ / # i n c I u de~< s t d d eTTi > 

//define stmpi(a) strupr(a) 
//define strlvvrfai 3 tr ' wr < a ) 
//define inp(ci} _inp.'a) 
//define inpw(a) snpw(a) 
//define outp(a.b7 outp(a.b) 
//define outpwfa bT _outpw(a b) 

//ifdef DO_WIN32 
//ifdef DO AD12 18 1 
//define DEBLEVEL 1 
//define DEBUG 

//include <debug h.> /*- DDK Debug Header Tile -7 

////define ■ Dehug_Printf Service LCODE DebugPrintf Service 

////include -'wdwiaps h> 
//endif I' DO_ADI2 181 7 
//define far 
//define near 

//define NULL ((void ~)0) 
//endif /' DO_WlN32 7 . 

//else /* DO_WIN 7 

tfifndef DO_CX018 A 
/'//include < null b^l 
flifndef NULL 

M definedf TINY ) || definedf SMALL_ ) || definedf MEDIUM ) 

//define NUEE 0 — 
//else 

#define NULL OL 

//endif 

//endif 

//endif /' DO_CX018_A 7 
//endif /* DO WIN 7 
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/ * * 

typedef unsigned char uin!8: /* range 0-255 7 

typedcf unsigned char umod8: /* 8 bit math requirements 7 

typedef unsiqned char octet. /* octet data 7 

#ifdef DOFORCED BYTE ALIGNED 

typedef unsiqned short oack*T 

#else /* OOjFORCED_BYTE_ALIGNED 7 

typedef unsigned char pack8: /" packed octet data */ 

ffendif /* DO^ FORCED BYTE^ALlGNED 7 

typedef packB expand8~7* expand cctet data to uintt6 7 • 

typedef signed char smt / . 

typedef unsigned short uint16: 
typedef signed short sint 1 5. 
typedef signed short sfractlS. 

typedef unsigned tony uint32. 
typedef signed lony sint?i 

ffifdef DO_CX025 
#define far 
#define interrupt 
#endif /* DO CX025 7 



/* 



#defme MAX UINT8 (255) 
#define MAX SINT7 (\27) 
^define MINJSINT7 (128) 

#define MAX UIN r 16 (65535) 
tfdefine MAX"~SINT15 (32767) 
#define MJNJ5INT15 (-32768) 

#define MAX UINT32 (4294967296) 
#define MAX"SINT31 (2 147483647) 
tfdefine MIN_T>INT3 t (-2147483648) 

tfdefmeMIN DOUBLE ( 1 7e-308) 
^define MA)C_UOUBLE (1 7e308) 



7 



typedef uint8 byte; 
typedef uint 16 word: 

typedef int bool. 

tfifndefDO CX025 
#define TAESE 0 
#define TRUE i.TALSE) 
#endif /' DO_CX025 7 

#if defined(DO ADI218D || defined{DO_ISARl 

typedef uint32 9sp_pm_t. /* Of size of DSP program memory contents 7 

typedef uint 16 dsp_dm t; /* Of size of DSP data memory contents 7 

typedef uint 16 dsp addrj: r Of size of DSP addresses 

typedef uint 16 dsp" xddr t. /* Of size of extended DSP addresses 7 
#endif /* DO ADI2T8 1 II DO ISAR 7 



#define forever for ( .) 
#ifdef DO CX025 

#define inrerrupt_disable_6nio stack() \ 

asm( 'ori w #0x0700. sr") 

#define interrupt_restore_from stack() \- 

asm( "andi w #0xF8FF. sr'*) 

#else/* DO__CX025 7 

#ifdefDO CX018 _A 

#define inferrupt disable onto stackQ \ 

{asm(" ORI #$700,SR\n NORn NOP\n").} 
tfdefine interrupt restore from stackQ \ 

(asm(" ANDI #$f8ff.SR\h NOP\n NOP\n").} 
#else/* DO_CX018_A7 

#ifdefDO WIN 

#define inrerrupt_disable_onto_stack() 



7 
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//define interrupt restore .from stack() 
#else r DO_.WIKT 7 ~ ~ 1 

Refine interrupt_disable_onto _slackM (asm pushf:asm cli;} 
//define interrupt restore from" stackf) {asm popf;} 
ttendit r DOym 7 

#endif r DO CX0I8 A 7 
#endif r DO3CX025" T 7 

#ifdef.DO CAPI20 

//define DEBLEVEL 1 

#pragma code_seg("_LTEXT'\ "LCUOE") 

include <debug.h> 

//include <vxdwraps h> 

#e(se /* DO CAPI20 V 
#ifdefDO I5AR 

typedef const chat far* LPCSTR. 
typedef char far' LPSTR: 
^define PASCAL pascal 
ffendif /* DO ISAR 7 
ftendif r DO* OAPI20 7 

' / 

#;fdef DO_WIN32 
////define sprintf Sprintf 
ttenrXif r DO WIN32 7 

/* * * , **** v 

/ v 



//define print_debucj off(A) 
define print diag oTf(A) 

//ifndef DO_CAPI20 

//define pnntpriniary(A) print routine A 

//define print_std(A) print routine A 

//define prinMnfo(A) prinf^routine A 

//define print_diag_6n(A> print_routine A 

//define print debug_on(A) print routine A 

//define report anomaly(A) report anomaly routine (A) 

//else /• DO_CAPI20 V " " 

//define print_std(A) Debug_Printf A 

//define print_;nfo{A) Debuq_ Printf A 

//define print_diag_on(Ai L^ebug^Printf A 

f/defme print debug_on(A) Debuq_Printf A 

//define reporr_anomaly( A) repOif_anoma\y_rouhne (A) 

tfendif /* DO CAP 120 7 

flifdef DO_J>RINT S_DISAULfcU 

//undef pnn\ primary 

//undef pnnCstd 

//undef print_ info 

//undef print debug 

//undef repon_anoma!y 

//define pfint primary(A) 

//define prinf_stcl( A) 

//define prin|_info{A) 

//define pfint debug(A) 

//define report anomafy(A) 

//endif /* DO_ PRINT S_ DISABLED 7 

#ifde?DO CX025 
#undef prfnt_primarv(A) 
//undef prin!_std(A) ' 
//undef printjnfo(A) 
//undef print debug(A) 
//undef reporl anomafy(A) 
//endif r DO_CX025 7 

/ ; — / 

#endif /* SYSTEMSTANDARDH V 
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(54) Method for determining the attenuation of a PCM signal over a digital channel 



(57) A method of determining digital channel atten- 
uation; comprising the steps of: receiving a known train- 
ing sequence of PCM codes, which PCM codes are 
subjected to the attenuation within the digital channel; 
quantizing the received known training sequence of 
PCM codes according to a predetermined thresholding 
procedure; identifying identical PCM codes created as a 
result of the thresholding procedure; and, determining 
the attenuation of the digital channel based upon the 
identification of identical PCM codes. A method is also 
disclosed for determining a digital channel PCM code 
transformation comprising receiving a known training 
sequence of PCM codes, which PCM codes are sub- 
jected to the PCM code transformation within the digital 
channel, quantizing the received known training 
sequence of PCM codes according to a predetermined 
thresholding procedure, and determining the transfor- 
mation of transmitted codes to those received . A 
method is also disclosed for improved echo cancellation 
in a communications network having an analog and a 
digital modem, comprising saving codes transmitted 
from the digital modem to the analog modem for echo 
cancellation, transforming, by a mapping table, codes 
transmitted from said digital modem to codes received 
by the analog modem, and, using the received codes as 
a reference signal for cancellation of echo. A method of 
improved spectral shaping using a transmit shaping 
transfer function in a communications network having 



an analog and a digital modem, comprising, transform- 
ing, by a mapping table, codes transmitted from the dig- 
ital modem to codes received by the analog modem, 
using the received codes for transformation to their lin- 
ear value equivalent representations, and, applying the 
linear value representations to the transmit shaping 
transfer function. 
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1. Claims: 1-3 

Method of determining digital channel attenuation. 

2. Claim : 4 

Method of echo cancellation. 

3. Claim : 5 

Method of spectral shaping. 
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